A new Fermi smearing approach for scattering of multi-GeV electrons by nuclei 
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The cross section for electron scattering by nuclei at high momentum transfers is calculated 
within the Fermi smearing approximation (FSA), where binding effects on the struck nucleon are 
introduced via the relativistic Hartree approximation (RHA). The model naturally preserves current 
conservation, since the response tensor for an off-shell nucleon conserves the same form that for a 
free one but with an effective mass. Different parameterizations for the inelastic nucleon structure 
function, are analyzed. The smearing at the Fermi surface is introduced through a momentum 
distribution obtained from a perturbative nuclear matter calculation. Recent CEBAF data on 
inclusive scattering of 4.05 GeV electrons on 56 Fe are well reproduced for all measured geometries 
for the first time, as is evident from the comparison with previous calculations. 
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I. INTRODUCTION 

Electron scattering by nuclei at high momentum trans- 
fers is a powerful tool to study the effective constituents 
of hadronic matter and exhibits very interesting new fea- 
tures. Within the Q 2 = -q 2 > l(GeV/c) 2 (or q = |q| > 
lGeV/c) domain, where q = (o>, q) is the four-momentum 
transferred by the virtual photon, the struck nucleon is 
relativistic having a momenta of the order its mass M. In 
addition for such regime the probability for exciting in- 
ternal degrees of freedom of the nucleon (nucleon inelas- 
tic response) becomes increasingly important. Ideally, 
to describe the target response, one should start from a 
relativistic covariant theory of nuclei. However, such an 
approach is not practicable due to the difficulties in treat- 
ing the meson exchange interactions. On the other hand, 
how to describe a nucleon with three momentum p inside 
the nucleus is a well known nonrelativistic nuclear struc- 
ture problem. Thus, a model that couples both regimes is 
necessary. Electron scattering experiments have been de- 
scribed with a great variety of approximations, starting 
with the plane wave impulse approximation PWIA. In 
the Born approximation the A(e, e')A' differential cross 
section reads 
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being k, k' = |k|, |k'|, W^ v the nuclear response tensor, 
L^(k, k') = l/2[k'»k v + k»k' v + (q 2 /2 - m 2 /2)g^} the 
lcpton tensor describing incoming and outgoing plane- 
wave electron states of four-momentum k = (e = 



Vk' 2 + m 2 , k) and k' = (e' = Vk' 2 + m 2 , k') respec- 
tively, and Q' = (9, 4>) the scattering angle. The PWIA 
lies on the following assumptions: 

i) the nuclear current operator can be written as the 
sum of the one-body nucleon currents; 

ii) the target decays virtually into a on-shell (A-l) nu- 
cleus (spectator) and the off-shell (p 2 ^ M 2 ) struck nu- 
cleon, of four-momentum p — (po, p) ; and 

iii) the nucleon that absorbes the photon is the same 
that leaves the target without interaction with the spec- 
tator, the final state interactions (FSI) being dropped. 
Under these suppositions, the nuclear response can be 
expressed as a convolution [1] 

q) = E / dEd P P m < (E, p)w% (p, g), (2) 

of the nucleon response w^{p,q) (m t = 1/2 and —1/2 
for protons and neutrons respectively) with the nuclear 
spectral function P ,nt (E,p). This gives the joint prob- 
ability of finding a nucleon with three momentum p in- 
side the target nucleus, and remove it with an energy 
E = E B + Ef^ v E B = Ma-i +M- M a is the nucleon 
binding energy and E A *^ 1 the excitation energy in which 
the residual nucleus is left. Notice that for an off-shell 
nucleon, the energy pq = po(E, p) depends on its remov- 
ing energy and its three-momentum, thus to implement 
the PWIA or any extension including FSI one must ad- 
dress some important questions. First, the nucleon struc- 
ture function is determined experimentally from proton 
or deuteron scattering on on-shell (free) targets, being 
p 2 = M 2 (or pq = E p = \Jp 2 + M 2 ). In our case we 
treat with an off-shell bounded nucleon with po ^ E p , 
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and po = po(E,p) depends on how the binding effects 
are included. Second, we need to extend the on-shcll mi- 
cleon structure function to the off-shell regime to use it 
as input in the nuclear response calculation. The min- 
imal hypothesis adopted in majority of works is to as- 

. n m+ (off-shell) / \ m+ (on-shcll) /~ , 

sume that w^ v [p, q) = w^v (p, q) , where 

p and q depend on the off-shell prescription adopted for 
Po = Po(E, p) • Third, whatever is the (p, q) pair we have 
a lack of the electromagnetic gauge invariance because 
^mt (off-shell)^,, ^ ^ o ^ e on _ gne jj ^ ff_sh e ii exten- 
sion. This brings in additional complications, a proce- 
dure being required to restore current conservation [1, 2]. 

The old data coming from the NE3 SLAC experiment 
[3] were analyzed within different approaches. One of 
the first PWIA calculations included onc-holc (lh) and 
two-particle - one-hole (2p-lh) excitations in P mt (E, p) 
[1, 4] . The cross sections, when expressed in terms of 
the well known x = Q 2 /(2Mlo) Bjorken variable, have 
been fairly well reproduced in the quasielastic peak [27] 
(x ~ 1) and inelastic (x < 1) regions, but underesti- 
mated for x > 1. In order to correct this the FSI were 
introduced in different ways . For instance, when the 
PWIA has been extended by assuming a factorization 
hypotesis for the final nucleus wave function and by in- 
troducing pair correlations [5], the discrepancies in the 
region 1 < x < 2 where circumvented. For x > 2 more 
than two nucleons should be involved in the scattering 
process and thus the use of an optical potential was re- 
quired [6] . Bcnhar et al. [7] have improved the PWIA re- 
sults by introducing the FSI through an optical potential 
and by generating a folding function from the multiple- 
scattering Glauber theory and color transparency. 
The quite recent CEBAF 4.05GeV electron scattering 
[8] experiment covers the range 1 < Q 2 < 7(GeV/c) 2 



and 0.2 



< „ < 



4.2, and vastly extend the angular and 



energy-loss range of the older NE3 SLAC one. Rinat and 
Taragin [9, 10] analyzed these results adopting an alter- 
native approach to the PWIA. The nuclear response func- 
tion is treated in a relativistic extension of the Gersch- 
Rodrigucz-Smith series [11], while the FSI were intro- 
duced through binary collisions. The CEBAF data are 
well reproduced for x < 1 and in the left hand side neigh- 
borhood of the quasielastic peak x> 1 [12]. Nevertheless, 
for all 9 geometries the calculated cross section overesti- 
mates the data by a factor up to 2-10 in the low energy 
lost region (x > 1), being these discrepancies associated 
to defects in the adopted momentum distribution. In 
the present work we develop a modified version of the 
PWIA (in the sense that the nucleon behaves as free one 
but with an effective mass) where the off-shell effects and 
FSI are included via the RHA, being at the same time 
the gauge invariance preserved. The Fermi smearing ef- 
fects are incorporated through a new nucleon momentum 
distribution, obtained from a perturbative calculation in 
nuclear matter. In addition, different parameterizations 
for the inelastic nucleon response measured at SLAC, are 
analyzed. The CEBAF data are satisfactorily reproduced 



for all measured geometries, taking into account that the 
cross section varies over many orders of magnitude, the 
mentioned overestimation being avoided in the x > 1 
region. 



II. ELASTIC AND INELASTIC CROSS 
SECTIONS 

As the electron probes a region of dimensions 1/q, 
for high momentum transfers, surface effects are sup- 
posed to be of minor importance and the nuclear mat- 
ter framework is adopted. How good is this assump- 
tion will be analyzed in Section IV, where the theo- 
retical results will be compared with the experimental 
and nuclear matter extrapolated data. The nucleon will 
be bounded by interaction with the scalar 4> and vec- 
tor mesons fields, within the framework of quantum 
hadrodynamics(QHDI)[13, 14]. The nucleus response 
tensor is the Lorentz invariant amplitude and reads [15] 
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(p'm' s m' u PLi\J(0UPA) 
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being P A = (M A ,0) and P{_ x = {yjpj + (Mj{_ 1 ) 2 , p/) 
the target and residual nucleus four-momentum respec- 
tively, with the mass M^_ 1 — M A -\+uj a _ 1 including the 

excitation energy u A _ v The sum on / encloses the set 
of final residual nucleus states. We also sum on the final 
states of the struck nucleon with four-momentum p' = 
(pg,p'), spin m' s and isospin m' t , with density V/(2n) 3 in 
the quantization volume V. J(x) is the effective hadron 
current density operator J^{x) = iip^T^ (q)ip(x) with 
T^q) = Fi{q 2 )l^ + iF2(q 2 )jM<J IJ , l ,q l/ for the nucleon elas- 
tic response case, being ^(x) and k the nucleon field and 
anomalous magnetic moment, respectively. 

We are going develop on the same footing the nu- 
clear response calculation within the mean field theory 
(MFT) (where the meson fields are approximated by 
their vacuum spectation, i.e. constant, values), and in 
the RHA[13] (where vacuum fluctuation corrections are 
added to the MFT results). Later, when we compare the 
calculated cross section with the data, the RHA election 
will be justified. The nucleon field is expanded as 
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where the single particle spectrum is given by 



(5) 



with E* = v/p 2 + M* 2 and M* = M + £(C S ,M*). 
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with p* = (E*,p) and [28] v* = p* ■ q/M* . n m *(p) is 
the nucleon momentum distribution in the target ground 
state |0a) 



n mt ( P ) 



V 
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pm,m t a P»i J mt 



,\0a), (8) 



normalized as 2 / dp n mt (p) = N m * , with N mt = Z,N 
for m t = 1/2, —1/2. The elastic Lorentz scalar functions 
present in (7) are 



FIG. 1: (a) Tadpole diagram included in the MFT and RHA 
self-energies, (b) Tadpole exchange diagram that is added 
in order to get the relativistic Hartree Fock self-energy. The 
dashed lines indicate the propagator of the scalar (S) or vector 
meson (V) that interacts with a nucleon n (full lines) 

M* < M is the effective mass acquired by the nu- 
cleon by action of the attractive scalar field and is deter- 
mined self-consistently [13] through the scalar self-energy 
£ = Smft ° r ^rha- ^mft includes the tadpole dia- 
gram a) in the Figure 1, retaining in its evaluation only 
the contribution from nucleons in the filled Fermi sea in 
the nucleon propagator (tick full lines). T,rha includes 
the same diagram but the full nucleon propagator (which 
encloses the contribution of the occupied negative-energy 
states) is used in the evaluation of the self-energy. Then 
the MFT or the RHA are derived by summing up the 
self-energy to all orders through the self-consistent de- 
termination of M* , being this procedure convergent in 
both cases. The first term in Eq.(5) accounts for the 
action of the repulsive vector field. Cy and Cs are the 
two free parameters [16], which depend on the meson cou- 
pling constants and masses, fixed to reproduce the exper- 
imental binding energy per nucleon of —16 MeV at the 
Fermi momentum pp = 1.42/m -1 (or the baryon density 
Pb = 0.19/to~ 3 ) for the normal nuclear matter. Assum- 
ing that the residual nucleus is left in its ground state and 
adopting the prescriptions ii) and iii) mentioned above, 
the response tensor can be obtained from Eqs. (3) to (5) 
as 



/M* 
dp—n m *(p)w™J(p*,q), (6) 
p 

where the factor 2 resembles the sum over spin states, 
and 
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where 
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F^(Q 2 )K r ' 



2M* 

(10) 
%t T and 
are the elec- 
tric and magnetic form factors, and r = Q 2 /4M* 2 . 
In the numerical calculations we adopt the Sachs 
form for them, assuming that they do not change 
in the nuclear medium [17]. Equations (7), (9) and 
(10) show that the MFT or RHA lead to the pre- 
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e 1,2 
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mt (on-shell) 



w. 



(p* , q) and 
(Q 2 ,^*), for the elastic 



case. The nucleon spinors carry a four momentum p* be- 
ing p* 2 = M* 2 , and as M* < M this makes us remember 
that the struck nucleon is bounded. Lorentz , parity and 
gauge invariances are now also fulfilled as were for a nu- 
cleon of mass M, as consequence of the form of the Eq.(7) 
[18]. FSI are included, since the nucleon is bounded 
also after the interaction with the photon. For Q 2 > 1 
(GeV/c) 2 the probability of exciting internal states of the 



nucleon is important, and a replacement w™\ 2 — > w 



1,2 



w 



+ w™* 2 m (7) should be done, adding an inelastic 
For w™* 2 we use different parametric 
fits done at SLAC for p(e, e')p' and d(e, e')d! data through 
the Eqs. (7), with M* = M. We assume that the recipe 



el,2 

contribution w™* 2 - 



o f f — shell 



(Q 2 ,v) 



-shell 



J il,2 



(Q 2 ,v*), which naturally 



appeares in the elastic case, is also valid for the inelas- 
tic nucleon response function. Finally, the decomposition 
w ei2 + w ii 2 leads also to split the inclusive cross 



w 



1,2 



section (1) in elastic and inelastic contributions. 



III. NUCLEON MOMENTUM DISTRIBUTION 

From Eqs. (2) and (6) it is clear that we are working 
within the FSA P m *(E,p) - n mt {p)5{E - E B ), giving 
n" H (p) the probability of finding a nucleon with momen- 



tum p, and isospin mt in the target |0a)- 



; (P) 



is 
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calculated in a OpOh + 2p2h + 4p4h configuration space 
for the A-target, being 



\0a) = A/" 



p's,h's 



+ 77ri2 C PlP2P3P4h 1 h 2 h 3 h 4 \PlP2P3P4hlh 2 h 3 h4) 
^ '> p'sJi's 

(11) 

where these \npnh), (with n = 0,2,4) stand for the un- 
perturbed states. In this way in the residual nucleus we 
have lh, 2p3h, 4p5h, lp2h,and 3p4h excitations when 
the struck nucleon is removed. The residual nucleon- 
nucleon interaction is included within a perturbative ap- 
proach as in Ref.[19] by expanding the coefficients C2 P 2h 
and C4 P 4h up to the first and second order, respectively. 
This "minimum" perturbative scheme allows to include 
norm corrections N = (0a|0a} _1 , avoiding in this way 
contributions of unbalanced disconnected diagrams. Wc 
get 



! (P) 



3N r 



47T P f, 



#(l-p) + 6n (2) (p)+<5n (4C) (p) 



(12) 



where p = |p| is measured in units of the Fermi momen- 
tum pp. The first term is the usual OpOh Fermi step 
function, while Sn^ 2 \p) and Sn^ 4C ^(p) (where the super- 
script C indicates "connected" 4p4h diagrams) enclose 
2p2h and 4p4h contributions respectively, which deplete 
it. The expressions for 5n^ 2 -*(p) and Sn^ 4C \p), are given 
in Ref. [19], while the Goldstone diagrams corresponding 
to them are shown in Figure 2. 



IV. RESULTS AND CONCLUSIONS 

We now compare the differential cross section calcu- 
lated within our model, with the CEBAF experimental 
results for 56 Fe [8] for the various accessible geometries 
6 = 15,23,30,37,45,55,74°. One of the parameteriza- 
tions for tu™* 2 was found by Bodek et al. [4] in the kine- 
matical range 1 < Q 2 < 20 (GeV/c) 2 and 0.1 < x < 0.77. 
The other one was reported by Whitlow [20], and cor- 
responds to the range 0.6 < Q 2 < 30 (GeV/c) 2 and 
0.06 < x < 0.9. 

The functions obtained in these parameterizations are 
described in detail in Ref. [14], and as they do not cover 
all the low energy and momentum transfer region of CE- 
BAF, an extrapolation is necessary. This fact could in- 
troduce some uncertainties in the calculation. Within 
the MFT and for 56 Fe (p F = 1.36 fin" 1 ), M* = 0.648M. 
This value is too low to reproduce satisfactorily the to- 
tal cross section since the quasielatic peak is shifted too 
much to the right and its width (Atv qe ) is enlarged in 



(a) 




(b) 



(c) 
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FIG. 2: Goldstone diagrams corresponding to the second- 
order 2p2h correction 5n^ (p) (a) and fourth order 4p4h cor- 
rection <5n' 4C ^(p) (b,c,d,e). Each line indicates schematically 
a particle or a hole state, the dots represent the residual in- 
teraction and the encircled dots correspond to the number 
operator n(p). With b, c, d and e we indicate different ways 
to attach the number operator to a particle or hole line in 
5n< 4C >(p). 



excess, as shown in the Figure 3. This behavior of the 
MFT at high momentum transfers was analyzed in Ref. [2] 
for the longitudinal response, where a dependence of <j> 
and V,j, with p was introduced [29] . The value of uj qe is 
certainly corrected (but not the detailed peak's shape) 
since as it is shown, the binding-energy shift effects are 
diminished as p increases. Nevertheless the p dependence 
of the fields carries a gauge invariance violation that is 
corrected by introducing a (nonunique) vertex correction 



The Eqs.(7), (9) and 

,,of J '-shell f n 2 



(10) 



are altered, and the 
prescription w" e [y°'™ l {Q 2 ,v) = w° e ^ shel1 (Q 2 , v*) is no 
longer valid. This brings a problem at the moment of 
introducing the inelastic response, since the assumption 



o f f — shell 



(QV) 



slid! 



(Q , v*), is not justified. 



Alternatively, we try to improve the MFT description by 
adding the vacuum fluctuation corrections to Y,mft [13], 
and go to the RHA [22] where M* = 0.74M. The Eqs. 
(7), (9) and (10), are now still valid. As can be seen 
in Figure 3 the RHA the binding-energy shift is more 
moderated and the width is diminished, getting a better 
description for the total cross section. This improvement 
is not casual since as it is well known the RHA yields 
to the "best" single-particle spectrum in the sense that 
it minimizes the energy of the whole system. In addi- 
tion when the longitudinal response for electron scatter- 
ing at q = 0.55GeV and l.lAGeV transfers is analyzed, 
we get a value uj qe = 0.182GeV (Auj qe = 0.302GeV) 
and 0.615GeV (0.436Gey) respectively, in fully agree- 
ment with the results reported in Ref. [2]. This indicates 
that keeping p-independent fields and changing the value 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



Elastic + Inelastic 

• Exp. 30° 
M=M 

M=0.648 M 

M=0.75M 



0.5 1.0 1.5 2.0 2.5 
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FIG. 3: We show the sensibility with the effective mass M* 
of the quasielastic and inelastic contributions to the cross sec- 
tions per nucleon for Fe. Here the replacement e' = e + lu 
is done. In the panel (a) both cross sections are shown sep- 
arately for the values M* = 1, 0.64, and 0.74. Thin lines 
indicate elastic cross sections while thick lines indicate the 
inelastic one. In panel (b) the total elastic + inelastic cross 
section is shown for the different values of M* . Again, exper- 
imental results come from Ref . [8] . 

of M* in Eqs. (9) and (10), one can still improve the 
quasielastic peak position both at low and high momen- 
tum transfers at the same time. FSI are taken into ac- 
count in our model at the RHA level. Binding effects are 
present in the final state, since the nucleon still has mass 
M* after absorbing the photon. This simple form of in- 
troducing FSI has never been used previously to describe 
a multi-GeV electron experiment with the inclusion of the 
inelastic nucleon response, being only described in the 
past the quasielastic cross section at intermediate ener- 
gies in the MFT framework [23]. 

FSI affect directly the quasielastic response defined in 
Eqs. (7), (9) and (10), since the size and position of the 
quasielastic peak are controlled by r = Q 2 /AM* 2 (which 
scales quadratically with M*) and v* , respectively. In 
the inelastic nucleon response, FSI effects are included 
indirectly through the replacement v — > v* = p*.q/M* 
in the tu™' 2 on-shell functions. Our results for the total 
cross section are shown in Figure 4. As can be seen in, 



co[GeV] 

FIG. 4: Calculated differential cross section per nucleon for 
different geometries for 56 Fe. Experimental data come from 
Ref. [8] . Results are shown for both, the fitting of Bodek and 
Withlow of the inelastic nucleon response, and for a value 
M* = 0.74 corresponding to the RHA. 



the overall agreement is good for all angles 9, considering 
that that the cross section varies over several decades. At 
ui < uj qe (x > 1) Withlow's fit seems to be preferred to 
Bodek's, which is due possibly to differences in the ex- 
trapolation for the x > 1 range. For u> > u> qe (x < 1) 
the behavior is opposite. We see that the model tends to 
overestimate the x > 1 data, in the last two 9 values. The 
inelastic response dominates the cross section at these ge- 
ometries since Q 2 ^ 4(GeV/c) 2 , and this overestimation 
could be also as consequence of uncertainties in the ex- 
trapolation for x > 1. We conclude that to implement the 
generated nuclear matter momentum distribution within 
the RHA framework, is a consistent approach for treating 
electron scattering at these momentum transfers. This is 
further supported by the following observations: 

i)In the Figure 5 we compare the nuclear matter ex- 
trapolated cross section reported in Ref. [3] for e = AGeV 
and 9 = 30° in the shown energy lost region, with the 
experimental results of Ref. [8] and our calculations. As 
can be seen the extrapolated results for this geometry 
are larger than the experimental data, being the differ- 
ence small comparatively to the range of variation of the 
cross section (in the Figure 5 we show only an interval of 
the full energy lost region). This indicates that nuclear 
matter is a reasonable framework at these momentum 
transfers. Our nuclear matter results change from below 
the data to above the extrapolation, indicating that some 
improvements as would be a momentum dependence of 
the effective mass or higher-order contributions to the 
momentum distribution should be more deeply analyzed; 
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ous theoretical evaluations of the cross section at x > 1 
[12], is not present in our calculation. This can be seen in 




co[GeV] 

FIG. 5: Comparison of the extrapolated nuclear matter (NM) 
cross section per nucleon at e = A.GeV reported in Ref.[3], 
with the experimental data at t = 4.05GeV [8] and our cal- 
culations, for the indicated geometry 

ii) We are working within the Fermi smearing approach 
that in previous calculations with other nuclear matter 
momentum distributions [24, 25], lead to a large overes- 
timation of the cross section at low electron energy loss 
[1, 26]. As can be seen from Eq.(12) and Figure 2, we are 
including fourth-order corrections to n(p) in addition to 
the usual second-order contributions [24]. In Ref.[19], it 
has been shown that the second-order perturbation ap- 
proach overestimates the depletion of the Fermi surface 
and thus the high momentum tail of the momentum dis- 
tribution. This could be the reason of the mentioned 
overestimation present in previous nuclear matter calcu- 
lations. In the Figure 6 we show the momentum distribu- 
tion obtained from the Eq.(12) together with its second 
order approach, being Sn^ 4C ^(p) dropped. In the same 
figure we show the momentum distribution of Ref.[24] 
(parameterized in [6]), which was obtained within a sec- 
ond order perturbation approach over a set of unper- 
turbed variational wave functions. It is important to 
mention that FSI also are responsible for the behavior 
of the cross section in the low energy lost region. In 
our model FSI are taken into account by using an ef- 
fective nucleon mass, and different values of M* lead to 
different contributions to the cross section in the men- 
tioned region, as can be seen from Figure 3. The value 
of M* = 0.74M within the RHA, seems to introduce FSI 
consistently with the implemented momentum distribu- 
tion; 

iii) Finally the overestimation by a factor 2-10 in previ- 
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1x10 3 - 




1x1C P 

12 3 4 

p[fm 1 ] 

FIG. 6: Comparison of the nuclear matter momentum dis- 
tribution used in our calculations (full lines), with its second 
order approach (dotted lines) and the momentum distribution 
of Ref.[24] (dashed lines). 

Figure 7 where we compare the ratios of the theoretical 
to experimental cross section, for some selected geome- 
tries where the differences are appreciable, in our model 
and in that of Ref.[12]. 

In summary, to treat the scattering of GeV electrons 
by nuclei we have implemented a new Fermi smearing ap- 
proach. Binding effects and FSI are introduced through 
the nucleon effective mass within the RHA, that leads to 
better results than the plain MFT [14] . In the model, cur- 
rent conservation is preserved naturally without ad-hoc 
modifications in the structure functions. Fermi smearing 
effects are introduced through a new momentum distri- 
bution that accounts for 2p2h and 4p4h correlations in 
the target, generated via a perturbative approach in nu- 
clear matter. We get a reasonable overall description of 
the behavior of the measured cross section at CEBAF, 
for the scattering of 4.05 GeV electrons on 56 Fe. The 
agreement for all the accessible geometries, has been sig- 
nificantly improved in comparison with previous theoret- 
ical studies [12]. It could also suggest that within the 
model the FSI and Fermi smearing effects combine con- 
sistently, this could be more clearly established examin- 
ing the effect of changing M* on the contributions to the 
cross section coming from the hole and particle strength 
functions. This more detailed analysis and the scaling 
behavior of the model will be reported elsewhere [21]. 
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